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Abstract. We prove global well-posedness for the microscopic FENE model under 
a sharp boundary requirement. The well-posedness of the FENE model that consists 
of the incompressible Navier-Stokes equation and the Fokker-Planck equation has 
been studied intensively, mostly with the zero flux boundary condition. Recently 
it was illustrated by C. Liu and H. Liu [2008, SIAM J. Appl. Math., 68(5):1304- 
1315] that any preassigned boundary value of a weighted distribution will become 
redundant once the non-dimensional parameter b > 2. In this article, we show that 
for the well-posedness of the microscopic FENE model (6 > 2) the least boundary 
requirement is that the distribution near boundary needs to approach zero faster 
than the distance function. Under this condition, it is shown that there exists a 
unique weak solution in a weighted Sobolev space. Moreover, such a condition 
still ensures that the distribution is a probability density. The sharpness of this 
boundary requirement is shown by a construction of infinitely many solutions when 
the distribution approaches zero as fast as the distance function. 
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1. Introduction 

It is well-known that the following system coupling incompressible Navier-Stokes 
equation for the macroscopic velocity field v{t,x) and the Fokker-Planck equation 
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for the probability density function f(t, x, m) describes diluted solutions of polymeric 
liquids with noninteracting polymer chains 

(1.1) dtv + {vV)v + Vp = V-Tp + uAv, 

(1.2) V-v = 0, 

2 2kT 

(1.3) dtf + {v-V)f + Vm-i'^vmf) = -V„ ■ (V„^(m)/) + — A^/, 

where x G M" is the macroscopic Eulerian coordinate and m G is the microscopic 
molecular configuration variable. In this model, a polymer is idealized as an elastic 
dumbbell consisting of two beads joined by a spring that can be modeled by a vector 
m (see e.g [1]). In the Navier-Stokes equation fll.ip . p is hydrostatic pressure, z/ 
is the kinematic viscosity coefficient, and Tp is a tensor representing the polymer 
contribution to stress, 



Tp = \ J m ® Vm'^{m)fdm, 



where \Ef is the elastic spring potential and A is the polymer density constant. In 
the Fokker-Planck equation (11.31) . ( is the friction coefficient of the dumbbell beads, 
T is the absolute temperature, and k is the Boltzmann constant. Notice that the 
Fokker-Planck equation can be written as a stochastic differential equation (see 
One of the simplest model is the Hookean model in which the potential is given 

by 



2 
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where H is the elasticity constant. A more realistic model is the finite extensible 
nonlinear elasticity (FENE) model with 



|2 



(1.4) ^(m) = — —hgll--^], meB. 

Here B -8(0, ^/h) is the ball with center and radius Vh which denotes the 
maximum dumbbell extension. In this work we shall focus our attention on the 
potential (11.41) and the case b > 2, which is known to contain the parameter range 
of physical interest. We refer the reader to [TJ H] for a comprehensive survey of the 
physical background. 

In past years the well-posedness of the FENE model fll.ip - (ll.3p has been studied 
intensively in several aspects. For local well-posedness of strong solutions we refer the 
reader to [13] for the FENE model (in the setting where the Fokker-Planck equation 
is formulated by a stochastic differential equation) with 6 > 2 or sometime 6 > 6, 
[9] for a polynomial force and [28] for the FENE model with b > 76. For a prelim- 
inary study on some related coupled PDE systems, we refer to the earlier work 
(however, the FENE model was not addressed there). Moreover, the authors in 
proved global existence of smooth solutions near equilibrium under some restrictions 
on the potential; further developments were made in subsequent works [T^IIIZ]. More 
recently, N. Masmoudi [21] proved global existence for the FENE model fll.ip -I IT^l) 
for a class of potentials (11. 4p with b > assuming that the data is small, or the model 
is restricted to the co-rotational case in dimension two. 
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For results concerning the existence of weak solutions to coupled Navier-Stokes- 
Fokker-Planck systems and a detailed survey of related literature we refer to [U 1201 
l23l El |3]. For an earlier result on existence of weak solutions, we refer to [8] for the 
Fokker-Planck equation alone with > 4. On the other hand, the authors in [T^ . 
investigated the long-time behavior of both Hookean models and FENE models in 
several special flows in a bounded domain with suitable boundary conditions. 

The complexity with the FENE potential lies mainly with the singularity of the 
equation at the boundary. To overcome this difficulty, several transformations relating 
to the equilibrium solution have been introduced in literature. See, e.g. [SI El El El 
[T3] . A detailed discussion will be given in Section 2. In [21], C. Liu and H. Liu closely 
examined the necessity of Dirichlet boundary conditions for the microscopic FENE 
model. By the method of the Fichera function they were able to conclude that b = 2 
is a threshold in the sense that for 6 > 2 any preassigned boundary value of the ratio 
of the distribution and the equilibrium will become redundant, and for b < 2 that 
value has to be a priori given. For the microscopic FENE model, singularity in the 
potential requires at least the zero Dirichlet boundary condition 

(1.5) f\dB = 0. 

This is consistent with the result in [12], which states that the stochastic solution 
trajectory does not reach the boundary almost surely. 

The boundary issue for the underlying FENE model is fundamental, and our main 
quest in this paper is whether one can identify a sharp boundary requirement so that 
both existence and uniqueness of a global weak solution to the microscopic FENE 
model can be established, also the solution remains a probability density. The answer 
is positive, and we claim that / must satisfy the following boundary condition 

(1.6) fd"^ \qb = for almost all t > 0, 

where d '==' d{m, dB) denotes the distance function from m E B to the boundary dB. 
Our claim is supported by our main results: the global well-posedness for the Fokker- 
Planck equation stated in Theorem |2l the property of the solution as a probability 
density given in Proposition |3l and the sharpness of (11.61) stated in Proposition |H 

In this article, we focus on the underlying Fokker-Planck equation (II. 3p alone. 
Let v{t,x) be the velocity field governed by (11.11) and (II. 2p . We assume that this 
underlying velocity field is smooth, then a simplification can be made by considering 
the microscopic model (II. 3j) along a particle path defined as 

dtX{t, x) = v{t, X(t, x)), X{t = 0,x) = X. 

For each fixed x, the distribution function f{t,m]x) '==' f{t,X{t,x),m) solves 

~ 2 ~ 2kT 

dtf + ■ (Vvmf) = -V^ ■ (V™^^(m)/) + —A„J. 

By a suitable scaling ([21]), and denote / still by f{t,m) = f{t,'m]x), we arrive at 
the following equation 

1 fbm \ 1 

(1.7) dtf + W- (nrnf) = -V ■ (^— / j + -A/. 
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Here, p = b — |mp and K,{t) = Vf (t, X(t, x)) is a bounded matrix such that Tr(/t) = 0. 
We omit m from in fH-71) for notational convenience. In this paper we prove 
well-posedness of (11. 7p subject to some side conditions. The well-posedness of the full 
coupled system (II. ip - (1 1.31) is the subject of a forthcoming paper |22j . 

A weak solution of the Fokker-Planck equation (II. 7p with the initial condition 

(1.8) /(0,m) = /o(m), m G 5, 

and boundary requirement (II. 6p is defined in the following. 

Definition 1. We say f is a weak solution of (ll.7p . (11.81) . and (11.61) if the following 
conditions are satisfied: 

For an arbitrary subdomain B' of B such that B' G B and almost all t G (0,T), 

(1) f eL\0,T-H\B')) anddtf eL\0,T-H-\B')), 

(2) for any^^Cl{B), 



:i.9) 
(3) 



dm = 0, 



(1.10) /(0,m) = /o(m) znL\B'), 
(4) and for Br =' 5(0, r), 

(1.11) lim \\fd-'\9Br \\LHdBr) = 0. 

Note that (ll.lOp makes sense since / G C([0, T]; L^(i?')) implied by (1) above, and 
also fd~^ \qb^ is well defined in L'^{dBr) by the standard trace theorem. 
Regarding the weak solution defined above, several remarks are in order. 

• The reason for taking compactly supported functions as test functions in Defi- 
nition 1 is that we want to avoid any priori restriction to a particular weighted 
Sobolev space. It is this treatment that allows us to prove sharpness of bound- 
ary condition (II. lip . 

• Boundary condition (II. lip ensures that f{t, ■) G L^{B) for each t. Indeed, we 
can choose ro G (0, a/6) such that if r > tq 

\\fd~^\aBr\\L-i{dBr) < 1- 

Then 

\f\dm = [ \f\dm+ [ \f\dSdr 



B JBr^^ J TO JdBr 

rVb 

< Ci\\f\\L2^Bro) + C2 \\fd~^\dBr\\LHdBr)dr<oo. 



ro 



Boundary condition (11.1 ip or (II. 6p is, in its type, different from the zero flux 
boundary condition 



(1.12) ('ir + ^^-^-'^"^^)-R 



= 0, 

dB 
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which is known to preserve the conservation property, and has been used in 
many priori works. The relation of these two types of boundary conditions 
will be discussed in Section 2 as well. 
In order to establish an existence theorem, we now identify a subspace of H^{B) with 
an appropriate weight to incorporate boundary requirement (11.111) . For simplicity, 
we consider the case with trivial velocity field such that k = 0, then equation (11.71) 
becomes 



It follows from this conservative form that the only equilibrium solution /'^'^ must be 
a multiplier of p^f^, i.e. 

where Z is a normalization factor such that f^'^dm = 1. 
We are interested in the case 

(1.13) b>2. 

In such f^'^ satisfies boundary requirement (11.611 . Moreover 

r e H\/,{B). 

Here, H\^^{B) = {0 : 0,(9^^.0 G L\^^{B)] with 

L\/^{B) = i^(t) : j 02(5_ |^|2)-V2^^< ooj. 

Our main results are summarized in Theorem [21 Proposition [3] and Proposition H] 
below. 

Theorem 2. Assume (11.131) and n{t) G C[0,T] for a given T > 0. 

(1.14) /o(m) G L\/,{B), 

then there exists a unique solution f of (ll.7p . (11.81) . and (ll.6p in the sense of Defini- 

tionUl Moreover, 

(1.15) 

O^aX ll/(^'OllL2^/,(B) + ll/llL2(0,T;Hi,/,(iJ)) + PJI|L2(0,T;^^ < | /o 1 1 L%/,(B) • 

(a) For any 

/o(m) G LUB), 
there exists at most one solution f . 

Proof. The proof of (i) will be done in Section 3-5. In order to prove (ii) we assume 
that /i,/2 are two weak solutions of the problem with arbitrary initial data fo{m). 
Then /i — /2 solves (ll.7p with zero initial data which is in L^_i,/2- From (I1.15P in (i) 
it follows that /i = /a in ^^(0, T; i^l^/s)- □ 

We remark that the restriction on b in (I1.13P is essential to obtain the energy estimate 
dLlS]). 

The weak solution thus obtained is indeed a probability density. More precesely 
we have the following. 
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Proposition 3. Let f be a weak solution to (11.71) . (II. 8p . and (ll.6p defined in Defini- 
tion 1 subject to condition (I1.14p . Then, 

(1.16) / f{t,m)dm= / fo{m)dm, Wt > 0. 

J B Jb 

Furthermore if fo{m) > a.e. on B, then f{t, ■) > a.e. on B for all t > 0. 
This proposition will be proved in Section 2. 

The following proposition states that boundary condition (II. 6j) is sharp for the 
uniqueness of the weak solution. 

Proposition 4. Assume (11.130 and K{t) G C[0,T]. If boundary condition (11.60 fails, 
that is, 

fd-' \dB ^ 

is assumed, then the Fokker-Planck equation (ll.7p with fo{m) = has infinitely many 
solutions. 

In other words, Proposition H] implies that part [ii) in Theorem [2] would fail if 
boundary requirement (II. 6p were weaken so that near boundary the distribution ap- 
proaches zero not faster than the distance function. 

The justification of sharpness follows from the existence of a Cauchy-Dirichlet prob- 
lem for w defined by 

(1.17) "^^p"^ 

with g being a class of functions properly constructed. 

This article is organized as follows. In Section 2, we prove Proposition [3] and provide 
some preliminaries including: (1) several transformations used to handle the boundary 
difficulty, (2) equivalence of two weighted function spaces, and (3) the relation of 
our boundary condition to the natural flux boundary condition. In Section 3, we 
transform the Fokker-Planck equation to certain Cauchy-Dirichlet problem, named 
as VT-problem, and define a weak solution of VT-problem in a weighted Sobolev space. 
The well-posedness of the py-problem is shown in Section 4 by the Galerkin method 
and the Banach fixed point theorem. This leads to the well-posedness of the Fokker- 
Planck equation. Theorem |2l details of the proof are presented in Section 5. In 
Section 6, we construct non-trivial solutions for the Fokker-Planck equation described 
in Proposition |H 

2. Preliminaries 

2.1. Probability density. With the definition of our weak solution given in Defini- 
tion 1 we shall show that / has the usual properties of a probability density function 
(i.e. it is non-negative and has a unit integral over B for all t > if it is so initially) 
- this is to prove Proposition [31 

Given /o in L'^_^i^{B) and /o > a.e. on B, we define /o,; = r]i * fo E C°°{B) for 
/ > 1. Here r]i{m) = l'^rj{lm) denotes the usual scaled moUifier. We have 

lim II /o,/ - /o||l2,„(b) = 0, fo,i > 0. 
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Suppose that /; is the weak solution of (11. 7p . (11.81) . and (11.61) subject to initial con- 
dition //(0,m) = fo,i{m) G C°°{B). Then, for any T > and < t < T, 



(2.1) 



f{t,m) - fi{t,m)dm 



< Cmax ||/(t,-)-/Kt,-)IL^/,(B) 



< C\\fo- fo,l\\L^_^^^{B)- 



Hence for justification of the conservation of polymers, it suffices to prove that 



(2.2) 



B 



fi{t,m)dm= / fi{0,m)dm, Vt > 0. 



To do so, we take a test function (p^ G C^{B) converging to xb as e — )■ such that 

< Vb-e 
>Vb-e/2 



\ 1, 


m 


lo, 


m 



and 
(2.3) 



From (ll.gp and the fact that derivatives of are supported in := B^_^i2\B^_^, 
we have 



(2.4) 



dtfupedm 



B 



B^ 



■ bfim-Vipe 1 

flKm ■ Vipe — -V// ■ \/(pe 



dm. 



Applying the mean value theorem of the form 

/ gdm = - / gdS, for some r ^ {Vb — e, Vb — e/2), 
Jb^ 2 Jg^^ 

to the first term on the right of (12.41) together with (12.31) . we obtain 



fiKm ■ VipedS 



dBr 



<C [ \fl\dS<C\\fid-'h2^aB^y 
J dBr 



Similarly the second term on the right of (12. 4p is bounded by 

dS<C\\fid-'\\L^9Br)- 



1 


fl 


/dBr 


p 



It follows from (II. lip that the above two upper bounds converge to zero as e — 
Integration by parts in the last term in (12. 4p yields 



B'^ 



< 



e 
2 

< C 



\fiA(ps\dm 



\fN¥^e\dS 



95= 



dBr 



\fiA^,\dS+ / 



dB^ 



\fl 



dS, 



dBriJdB^ 



which, in virtue of \fi\/s < \fi\d ^ on OB^ U is converging to zero as e — )■ as 
well. 



HAILIANG LIU AND JAEMIN SHIN 



Due to Theorem [2] and the initial condition Jq i G C°°{B), it follows that d^fi is 
bounded in any Br ioi < r < y/b. Thus, for any r, s > 



fi{T,m)ip^dm- I fi{s,m)ip^dm 



B 



B 



" d 

dt 



fi{t,m)ip^dm ] dt 



B 



dtfi{t,m)ipsdmdt 

B 

Using the estimate for dtfnij, m)(pi.dm together with the boundedness of fi{t, m), 
we can send e to zero to obtain fl2.2p as claimed. 

We now turn to justify the positivity. Consider the transformation introduced in 

m 



(2.5) 
Then wi solves 



(2.6) 
where 



p^dtWi - ^p'^Awi - p[ ^'^^ ^ rn - pKm] ■ Wwi - c{m)wi = 0, 



c(m) = —Kp^ + a[nh + (2a + 2 — n — h)\mf'] + (6 — 2a)pm ■ Km. 



Then for any B.,. C B, Wi is a classical solution in (0,T]x B^.. It was shown in [21] that 
there exist a < 6/2 — 1 and K so that c{m) < 0. The maximum principle yields that 
Wi can not achieve a negative minimum at the interior points of [0, T] x Bj.. Thus the 
negative minimum of wi, if it exists, can only be attained on the parabolic boundary 
of the domain. 

From the transformation (12. 5p and the condition b/2 — a > 1, it follows that the 
negative minimum of //, if any, can only be attained at the initial time. Therefore 

fi> - max i > 0. 

Now fix t. For any xq and r] such that B{xo', r]) C B, 



fdm+ fidm< / \f-fi\dm< / < C| |/o-/o,/| 

B(x-o;r?) JB{xo;ri) J B{xo;n) JB 

Here (12. 2p has been used to obtain the last inequality. Hence 



6/2 



fdm > 



' B{xo;ri) 

which as / — 7- oo leads to 



fidm - C\\fo,i - /o||l2 



B(xo;r]) 



6/2 ' 



[ />o. 

JS(xo,r,) 



Since xo and t] are arbitrary, / > almost everywhere on i? for t > 0. The proof of 
Proposition [3] is now complete. 
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2.2. Transformations. To overcome the difficulty caused by the boundary singu- 
larity, several transformations have been introduced in literature. With boundary 
condition (11.61) . in this work we introduce 



to transform the Fokker-Planck equation to a degenerate parabolic equation with zero 
boundary condition (see details in Section 3). A widely accepted transformation is 
the ratio of the unknown to the equilibrium solution, i.e.. 



Such a transformation was used in [21] to reformulate the Fokker-Planck equation, 
and examine whether a Dirichlet type boundary condition is necessary. 
A third transformation is 



This was used in [8], [H] to remove the singularity at the boundary in the resulting 
equation. It was also used in [15] to formulate a weak formulation of w for discretiza- 
tion using a spectral Galerkin approximation. 
Another transformation defined by 



with b > 4s^/(2s — 1) and s > 1/2 is said to also lead to a well-posed problem. 
The minimum value of the function 4s^/(2s — 1) is attained at s = 1, yielding the 
maximum range of b values, 6 > 4. This transformation was proposed in [S], E] in 
the special case s = 2 and s = 2.5, where these values were chosen on the basis of 
numerical experiments in two and three dimensions, respectively. We note that our 
transformation w = f/p corresponds to s = 1, but not limited by 6 > 4. 

2.3. Weighted Sobolev spaces. In contrast to the standard weighted Sobolev space 
H\i^{B) used in this work, the following weighted function space 



is well known in literature for Fokker-Planck equations with FENE potentials, see 
e.g. [H HH El [20], EU [15]. We now show their equivalence as long as > 2. 

The key estimate we need to prove the equivalence is the embedding theorem stated 
in Lemma [5l Set if) = 0p~^/^. If G H\^^i^[B)^ we use the relation 






2m 



pb/2 



,b/2+l 



It is obvious that 



V0 
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Also the use of Lemma [5] and the fact that H^^i^{B) = H\^i^{B) for 6 > 2 (see [T6 



-6/2'' 



give 



06/2+1 



< C 



6/2 



-b/2-2 -i)/2 

Hence G p^l^Hlj^{B). If e p^/^Hlj^{B) we use the following identity 
It is easy to see that p^/^V"?/^ G 6 > 2 we have 

llp^''^VlU,,,, = II^IUg,,_,<c^ll^ll.,V2 

by Lemma O Thus G H^i^^B). These together verify that p'^^'^Hli^ and H^^^^^ are 
equivalent when b > 2. 



2.4. Boundary conditions. Granted certain smoothness of /, e.g. / G C^{B), one 
may argue that our boundary condition (11.61) is equivalent to the zero flux boundary 
condition ffLT2|l . 

Set = and g = fd~^. We calculate the flux 



J :-- 



KUlf + - V/ 

2p 2 

b\m\ f 1 df 



m 
\m\ 



\m\v ■ Kuf. 



2{\m\ + Vb)d 2du 

Due to singularity on boundary it is necessary that f\QB = 0. For any point p G dB, 
let m be a point in B such that m + du = p. Then 

/(m) 



df, . 



We thus have 



lim 

d~^0 d 

1 



Jip) = lim J(m) = -(b-2) lim 

d-s>0 4 d^O d 



For b ^ 2, this implies that J(p) = if and only if 

/(m) 



lim 



0. 



d-5.0 d 

3. Transformation of the microscopic FENE model 

In what follows we shall call the Fokker-Planck equation (11.71) with initial condition 
(II. 8p and boundary condition (II. 6p as the Fokker-Planck-FENE (FPF) problem. We 
first formulate a time evolution equation from the FPF problem. Define w{t,m) as 

(3.1) f{t,m)=w{t,m)p. 

Then (ll.7p is transformed to 



(3.2) 

where 
(3.3) 



dtw 



-Aw 



{b — 4)m — 2pKm 
2p 



■ Vw 



-w 



P 



0, 



c{t, m) = 2m ■ K,{t)m + n{b/2 — 1). 
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Setting a parameter 
we rewrite (13.21) as 

r . 1 . . .J 

0. 



dtwp^ — -V ■ (Vwp^) + Km ■ Vwp^ — cwp^ ^ 



The boundary condition (11. 6p implies that w(t, ■) satisfies a homogeneous boundary 
condition for almost all t since the distance function d and p are equivalent (see (13. 9p ). 
The FPF problem is formally transformed to the following VT-problem: 

(3.4) dtwp^ - ■ {Vwp^) + Km ■ Vwp^ - cwp'^-^ = 0, {t, m) e (0, T] x B, 

(3.5) w{0,m) = Wo{m), m E B, 

(3.6) w{t,m) = 0, {t,m) G [0,T] x 95. 

Here, 

Wo{m) = /o(m)p^^ 

according to the transformation (13. ip . 

In order to define a weak solution of TV-problem we introduce a weighted Sobolev 
space H^{Q] a) for a nonnegative measurable function a as a set of measurable func- 
tion (h such that 



m{n;a) ■= / (I + 0')arfm < CX). 



Similarly, a weighted L'^{Q;a) can be defined. H^{Q;a) denotes a completion of 

o 

C^{n) with II 

■ l|Hi(Q;cr)- It is obvious that H^{Q]a) and H^{Q;a) are Hilbert spaces 
with the inner product (-, ■)H^(n;a) defined as 



^'l, 4>2)m{n;a) = / (V01 ■ V02 + 0102)o"C?"^ 

Jq 



and 

H\n;a) C H\n;a). 

For notational convenient, we use i?^(fi), ^),(^^) and L^(n) for if^(n; p''), H^{n; p^) 
and L^(f2; p^) respectively. We also omit the domain VL if it is obvious. 

Lemma 5. Suppose that Q = B. 

(1) IfcpeHl for p< 1, then 

(3.7) ||0lk^_, <Co||0||Hi. 

If p > 1, we have the same inequality for cf) G H^^ 

(2) If (p E for /i < 1, then the trace map 

4> ^ 4>\dn 
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is well defined, i.e. it is a hounded linear map. 

o 

In particular, for cf) G H^^ 

(3.8) r(0) = 0. 
Proof. In [25] (see also [16]), it was proved that 

provided dQ is Lipschitz continuous. Recall that d denotes the distance from m to 
the boundary of Q. fl3.7p follows from 

(3.9) Vbd< p< 2Vbd. 

It is also known that the trace map T is well defined for < /i < 1 ([251 HB])- For 

m\m<b-^/'m\Hi^ 

since p'^ > for all m G 5. Therefore, T is well defined for /i < 1. (13. 8 p is obvious 

o 

from the definitions of the trace map and if^. □ 

We now define a weak solution to VT-problem in a standard manner. Multiplication 
by a test function G C^{B) to the equation (13.41) and integration over B yield 



dtW(fp^ H — Vw ■ V(fp^ + Km ■ Vw^p^ — cwipp^ ^ 



dm = 0. 



2 

This equation is well defined assuming that dtw(t, ■) G {H^)* , the dual space of if^, 

o 

and ty(t, ■), G H^p due to the boundedness of c and Lemma |5l Moreover, 

o 



implies 

o 



C c (M)*. 



Thus 

«;(t,x) gC([0,T];L2). 

Here we identify with its dual space. 

Let {■,-)h denote the paring of a Hilbert space H with its dual space H* and 

(3.10) 'h[w,(f]t\ = - I Vw(t,m) ■ Vipp^dm + f nm ■Vw{t,m)(pp^dm. 

2 J B Jb 

We now describe the weak solution we are looking for. 
Definition 6. A function w{t,m) such that 

w{t, m) G L^{0, T; ^J), with dtw{t, m) G ^^(0, T; (^J)*) 
is a weak solution of W -problem, (I3.4p - (l3.6p . provided 

o 

(1) For each ip G i?^ and almost every < t < T , 

{dtw{t, ■),ip) o ^ + L[w, ip:,t] = / cw{t,m)ipp^~^dm, 
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(2) w(0,m) = wo{m) in sense, i.e. 

\w{0, m) — WQ{m)\^ dm = 0. 



The following energy estimate for 'L[w,w.,t\ for fixed t is crucial. 

Lemma 7. There exist positive constants Ci and C2 depending only on h and | |k| |l°°(o,t) 
such that 

Ci||«;(t,-)||^i<i[^,^;t] + C2|Kt,-)lli.. 

Proof. Let (p = w in (13.101) and apply the Schwarz inequality we arrive at the above 
estimate as desired. □ 

The well-posedness of the W-problem is stated in the following 

Lemma 8. W -problem, fl3.4p -f lXUl) . is uniquely solvable in weak sense for Wq E L'^p. 
Furthermore, 

may: \\wit, ■)\\t2 + \\w\\ o + ||(9/w|| o < C||woll/-2. 

0<t<T" ^'^11-^/3 II llL2(o,T;/fi) II IIl2(0,T;(H1)*) - H '^H^/S 

A detailed proof will be presented in next section. 

4. WELL-POSEDNESS for THE TRANSFORMED PROBLEM 

In this section, we show the well-posedness of the weak solution to VT-problem. For 
this aim, we consider the following [/-problem containing a non-homogeneous term 
/i(t,m)eL2(0,r;L2_^). 

(4.1) dtupl^ - ^ V ■ (Vwp^) + Km ■ Vup<^ - h = 0, (t, m) e (0, T] x B, 

(4.2) u{0,m) = Uo{m), mEB, 

(4.3) u{t, m) = 0, {t, m) E [0, T] x dB. 
The weak solution of [/-problem is defined similarly. 

Definition 9. We say a function u such that 

u e L'^iO, T; Hi), with dtU G ^^(0, T; (^J)*) 
is a weak solution of U -problem provided 

o 

(1) for each ip G and almost every <t <T 

(9fM(t, ■), V9) o _^ -|- L[m, yj; t] = / h{t,m)ipdm, 

Jb 

(2) n(0,m) = UQ{m) in L^. 

We remark that / h{t,m)Lpdm is finite for any h{t, ■) G L2_^ since (p G L'^^_2 from 
Jb 

(13. 7p . Thus / h{t, m)Lpdm can be understood as the Lq inner product although h{t, ■) 
Jb 

. +^ r ^ 



may not belong to Lq. 



The well-posedness for [/-problem follows from the standard Galerkin method. 
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Lemma 10. For given h G L'^{0, T; L\_p) and Uq G L^, U -problem has a unique weak 

solution. Moreover, 

(4.4) 

Proof. We first construct an approximate solution in a finite-dimensional space. Let 

o 

{(pi} be a basis of if^ and L^. The existence of such a basis can be verified from the 

I 

fact that is a dense subset of L^. Consider an approximation ui{t, m) = d[(t)(l)i, 

i=l 

where satisfies 



(4.5) {dtUi{t,-),(j)j)°,+Hui,<Pj;t] = {h{t,-),(f),)L2, 1<J<1, 

I 

(4.6) ^ d\{0)(l)i -)■ uq in L^, as / oo. 



i=l 



Since (14. 5 p and (14. 6 p form a system of linear differential equations, {(i-} is uniquely 
determined for each /. We rewrite (14. 5 p as 

(4.7) {dtUi{t,-),(f)j)L2+L[ui,(l)j;t] = (/i(t,-),0i)L2, 1 < J < /. 

Apply (i^- to (14. 7p and sum for 1 < j < I, then for almost every t 

{dtuiit, ■), + L[u/, u;; t] = {hit, ■), ^^K^^ ■))Lg- 

From Lemma [71 it follows that 
d 
dt 



(4.8) -||nKt,0llk+2C'i||nKt,0llm <2C'2||nKt,OII^ +2(/i(t,-),uKt,-))L^ 



KMt, ■),«Ki, ■))l?J <;ttIIM^,-)IIL +oColh(^,-)llk- 



From (13. 7p . for any 5 > 

^IIMV)lli5_,4 

With 5 = Ci/Cg, fOj) can be rewritten as 

(4.9) |||«Ktr)|li.+C'i||nKt,OII^. <2C2||nKt,0lli.+C'oVC'i||Mtr)H^ 
or 

|lK(t,-)lli^<2C2iK(t,-)lli.+CoVciiiMt,-)lli|_^. 

Use Gronwall's inequality to obtain 

m^^\\ui{t,-)\\l2 < C (||Mo|li| + ll^lli2(o,T;L2^^)) : 

where C is an appropriate constant which depends on /3, 6, T and On the other 
hand, integration of (14. 9 p from to T together with above inequality yields 



(4.10) o <c{\\uq 



|2 



L2(o,r;Hi) V" ""-^^ IIl2{0,T;L2_^) 
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A similar argument to that in [TU] gives us the estimate for as 

o 

Here we have used (14.51) with G Hj^ such that < 1 and f l4.10p . By passing to 

the hmit as Z — )■ oo and a standard argument (e.g. see [10]), we have well-posedness 
for [/-problem. □ 



Now, we introduce a linear map A to connect W and U -problems as 

A: L\0,T-Ll) ^L\0,T-Ll_p) 

Since c is bounded, 

\\cw{t,-)p^-'\\l, < \\c\\l^ I w\t,-)p'^~'p'~^dm, 

2-/3 



B 



|2 



w'^{t, ■)p^dm. 



B 



Thus, A is well defined and 

ll^(Wl) -^(w2)|li2(o,T;Lt^) ^ -W^2|li2(o,r;L2)- 

We define another map J-" such that 

7: C([0,r];L2) ^C([0,r];L2) 
w y-^ u. 

Here, J^{w) is given by the weak solution of [/-problem with 

h = A{w), 

and the initial condition 

Uo{m) = w{0, m). 

The map J-" is well defined from Lemma [TO] and the definition of A. Now we show 
that J-" is a contraction mapping for sufficiently small r. Let 

From the energy estimate (14. 4p . 

11^1 -«2||c([0,r];L|) < C\\A{Wi) - A{w2)\\l2^^^^,.Ll^^) 



C / \\A{w,){t,-)-A{w2){t,-)\\i._dt 



< C / \\wi{tr)-W2{t,-)\\i2dt 
Jo 

= Cr||u;i -^^;2||c([o,r];L2)• 



Thus, J-" has a unique fixed point w in C([0,r];L^) and w solves W^-problem in a 
weak sense in (0, r] x B, if Cr < 1. We are able to continue this procedure to obtain 
the global well-posedness for the above constant C is independent of r. 
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For the fixed point w, fl4.4p and the boundedness of A imply that for t G [0, r'] 

max Ollz-a + llif II °, +||(9tty|| 

0<t<T'" ^'^11^,3 II llL2(o,r';Hl) II IIl2(0,t';(M)*) 



r 2 



<C\\A{w)\\L2^Qy.^_^) + C'I^oIIl^ 

- O^t <r' 1 1 ^ ' ■ ) 1 1 ^1 + 1 1 ^0 1 1 l2 • 

We select a small t' <T such that Cr' < 1. Then 

max ||w(t, ■)||r2 + ||w|| o + ||c?jw|| o <C||wol|r2. 

Thus, 

lk(r',-)llL2 < C||tyo||L2 

and 

max ||w(t, ON r2 + ||w|| + ||9tw|| o < C||w(r', ■)ll7'2 

r'<4<2r' " ^ ^ " -3 L'^{t' ,2t';HI) H ^ ^ L^t' ,2r' ^{HI)*) " H ^ ' ^H^/? 

< C'^lko||L2. 

Continuing, after finitely many steps we obtain an energy estimation similar to fl4.4p . 
The proof of Lemma [H] is thus complete. 



5. Well-posedness for the FPF problem 

In Section 2, we transformed the FPF problem to W-problem formally, but it is 
not difficult to show that they are equivalent. Indeed, one can verify that boundary 
condition (11. 6p in the sense of (11. lip for the FPF problem is equivalent to the null 
boundary condition for py-problem. 

For any test function ip & C^, the weak solution formulation for / can be trans- 
formed to the weak solution formulation for w, with Lfp^^"^'^ as the test function. This 

o 

is valid since (pp^^"^'^ G is dense in if^. Such a justification can be reversed, hence 
the FPF problem and VT-problem are equivalent. 

Now we seek the function space in which the weak solution / to the FPF problem 
belongs. Recall that f3 = -b/2 + 2. For fixed t G [0,T], (jSTD implies 

(5.1) / Ifl^p-'/^dm = [ \w\^p^dm, 

Jb Jb 

(5.2) /(I/P + iV/Hp-'/'dm < C [ {\w\^ + \Vw\^)p^dm. 



B JB 



Also, for (j) G H^^i^ we have 

(5.3) mf,ct>)H^_J = |(9iii;,p-V)//^l < C||9,i/;||(^i). 11011^1^^^. 



The estimate of the weak solution, (11.150 follows from Lemma IH] together with (15. ip - 
(15. 3p . This finishes the proof of (i) of Theorem [2l 
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6. NON-UNIQUENESS 

In this section we show that fll.6p is sharp in the sense that more solutions can be 
constructed if a weaker condition is imposed — this is to prove Proposition HI 

It suffices to construct more than one solution to the Fokker-Planck equation with 
/o(m) = and the assumption 

(6.1) Wfd-^ \oB, Wi^dBr) T^Oasr^Vb for t e /. 

Here J is a nonzero measurable set. The idea is to consider a class of functions 
g{t,m) e W^'°°{{0,T) x B) such that g{0,m) = and g{t,m)\aB 7^ for t > (e.g. 
g{t,m) = t|mp) and show that for each g the following problem has a solution. 

1 fbm \ 1 

(6.2) dtf + V- i^mf) = -V ■ i—f\ + - A/, in (0, T] x B, 

(6.3) /(0,m)=0, meB, 

(6.4) f{t,m)p~^ = g{t,m), m{0,T]xdB. 
Note that /3 = —6/2 + 2 < 1, we can choose a parameter 7 such that 

(6.5) max{/3, -1} < 7 < 1. 
To proceed, we define 

/ 

w = g. 

P 

The resulting equation when multiplied by p^~"' leads to the following 

(6.6) dtwp'^ — ■ (Vwp^) + (/3 — 7)m ■ Vwp'^~^ + nm ■ Vwp^ — ho = 0, 

(6.7) w{0,m) = 0, me B, 

(6.8) w{t,m) = 0, {t,m) G [0,T] x dB, 
where 

hoit, m) = cwp"'^^ — dtgp'^' + -V ■ iygp^) — {13 — j)m ■ Vgp^~^ — nm ■ Vgp^ + cgp"*"^ 
with 

Let 



c(t, m) = 2m ■ K,{t)m + n{b/2 — 1). 



A: i^'(0,r;L2) ^^^(Cr;^^,) 

This is well defined since 7 > — 1 from (16.51) and the assumption that g e W'^'°°{{0, T) x 
B). From the same argument as that in Section 4, it follows that f l6.6l) -( l6l8l) has a 
unique solution w such that 



w{t,m) e L^{0,T;H}^), dtw{t,m) G L\0,T; {H}^)*), 
provided the corresponding ^/-problem 

(6.9) dtup"^ - IV ■ [Vup^) + (/3 - 7)m ■ Vup^'^ + Km ■ Vup^ - ho 

(6.10) u{0,m) = 0, m G 5, 

(6.11) M(t,m) = 0, (t,m) G [0,T] x dB, 
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has a solution for any G L^(0,T; L2_^). Note that 7 < 1 is essential in order that 
the trace of w at the boundary is defined. Equation (16. 9p is of the form of (14.11) but 
with an additional term (/3 — 7)m ■ Vup'^~^. We thus define 

Lo[m, (p; t] - I Vm ■ Vipp^dm + (/3 — 7) j m ■ Vu(pp^~^dm + f Km ■ Vu^pp^ dm. 

2 J B J B J B 

We may obtain the existence and uniqueness for (I6.9I) - (I6.11I) from the same argument 
of the well-posedness for [/-problem (I4.ip - (l4.3p . if there is an energy estimate of 

o 

Lo[M,M;t] which is similar to L[m, in Lemma [71 Indeed, for u G L^{0,T; H}^) 

- J \Vu\'^p'^dm = Lo[M,M;t] — ^ ^ ^ m ■ Vv^ p'^^^dm^ — j Km ■ Vuup'^dm. 
We now claim that 

(6.12) / m ■ Vu^p^-^dm < 0. 

J B 

Given this together with 7 > /3 from (16. 5p we have 

- I \Vu{-,t)^ p^ dm < Lo[u,n;t]— j nm-'Vuup^dm 

< Lo[n, n; t] + ||k||loo(o,t)V^ j \V p'^ dm + ^ j \u\^p^d'n 

for any 5 > 0. By taking 5 > \ \K,\\L°°{o,T)Vb, we obtain 

C;||n(t,-)||^i <Lo[u,u;t]+C^||n(t,-)||i2 

for appropriate constants C[ and Cg. 

To verify the claim (I6.12p . we define the trace operator To such that 



To: H}^{B) ^L\dB) 

U ^ up 2 \qb- 



Integration by parts on (16.120 yields 

I m ■ Vu'^p^-^dm = - I {np^-^ + 2(1 - 7)|m| V"^) dm + Vb [ u^p^-^dS, 

J B J B JdB 



or 



Vb / u^p^-^dS = 2 m- Vuup^'^dm + I {np^''^ + 2(1 - 7)|m| V"^) dm 

JdB Jb Jb 

< C\\u\\l,^. 

o 

Thus To is well defined, and for u G Hl^, %{u) = 0. Finally we obtain 

/ m ■ Vu^p"<~^dm = - I u'^ {np"<~^ + 2(1 - 7)|m| V"^) dm < 0. 
Jb Jb 

This shows that there is a unique weak solution u of (l6.9p -( l6rTTl) . and thus w G 

L\0,T;H]) of (S-S- 

Finally, f = {w + g)p is a solution of (I6.2p - (l6.4p satisfying (16. ip for each g. Hence 
the uniqueness of (I6.1l) - (l6.3p fails as stated in Proposition HI 
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7. Conclusions 

In this paper, we have identified a sharp Dirichlet-type boundary requirement to 
estabhsh global existence of weak solutions to the microscopic FENE model which 
is a component of bead-spring type Navier-Stokes-Fokker-Planck models for dilute 
polymeric fiuids. Such a boundary requirement states that the distribution near 
boundary approaches zero faster than the distance function. With this condition, 
we have been able to show the uniqueness of weak solutions in the weighted Sobolev 
space H^j^^^{B), which is equivalent to the widely adopted weighted function space 

p^/'^Hli^{B) for Fokker-Planck equations with the FENE potential. Moreover, this 
condition ensures that the distribution remains a probability density. The sharpness 
of the boundary condition was shown by construction of infinitely many solutions 
when the boundary requirement fails. In other words, such a condition provides a 
threshold on the boundary requirement: subject to this condition or any stronger 
ones incorporated through a weighted function space, the Fokker-Planck dynamics 
will select the physically relevant solution, which is a probability density, see e.g. 
[H HH El EHl EU [15], and converges to the equilibrium solution 

^-ipV2 [13]; any 

weaker boundary requirement may lead to more solutions, each depending on the 
rate of f/d near boundary. A detailed elaboration of boundary conditions for the 
coupled Navier-Stokes-Fokker-Planck model will be the goal of our work [22] ■ 

Acknowledgments 

Shin thanks Professor Paul Sacks for stimulating discussions on Proposition HI We 
thank the referee for valuable suggestions and pointing out two relevant references [15] 
and [3]. Liu's research was partially supported by the National Science Foundation 
under Kinetic FRG grant DMS07-57227 and grant DMS09-07963. 

References 

[1] J. W. Barrett, C. Schwab, and E. Siili. Existence of global weak solutions for some polymeric 

flow models. Math. Models Methods Appl. Sci., 15(6):939-983, 2005. 
[2] J. W. Barrett, and E. Siili. Existence of global weak solutions to kinetic models of dilute 

polymers. Multiscale Model. SimuL, 6:506-546, 2007. 
[3] J. W. Barrett, and E. Siili. Existence of global weak solutions to dumbbell models for dilute 

polymers with microscopic cut-off. Math. Mod. Meth. Appl. Sci. , 18:935-971, 2008. 
[4] R. B. Bird, C. Curtiss, R. C. Armstrong, and O. Hassagcr. Dynamics of Polymeric Liquids, 

Volume 2: Kinetic Theory. Wiley Interscience, New York, 1987. 
[5] C. Chauviere, and A. Lozinski. Simulation of complex viscoelastic flows using Fokker-Planck 

equation: 3D FENE model. J. Non-Newtonian Fluid Mech. , 122:201-214, 2004. 
[6] C. Chauviere, and A. Lozinski. Simulation of dilute polymer solutions using a Fokker-Planck 

equation. J. Comput. Fluids., 33:687-696, 2004. 
[7] M. Doi and S. F. Edwards. The Theory of Polymer Dynamics. Oxford University Press, Oxford, 

1986. 

[8] Q. Du, C. Liu, and P. Yu. FENE dumbbell model and its several linear and nonlinear closure 
approximations. Multiscale Model. SimuL, 4(3):709-731 (electronic), 2005. 

[9] W. E, T. Li, and P. Zhang. Wcll-posedncss for the dumbbell model of polymeric fluids. Comm. 
Math. Phys., 248(2):409-427, 2004. 
[10] L. C. Evans. Partial Differential Equations, volume 19 of Graduate Studies in Mathematics. 
American Mathematical Society, Providence, RI, 1998. 



20 



HAILIANG LIU AND JAEMIN SHIN 



[11] Y. Guo. The Vlasov-Maxwell-Boltzman system near Maxwellians. Invent. Math., 253:593-630, 
2003. 

[12] B. Jourdain and T. Lelievre. Mathematical analysis of a stochastic differential equation arising 

in the micro-macro modelling of polymeric fluids. In Probabilistic methods in fluids, pages 205- 

223. World Sci. Publ., River Edge, NJ, 2003. 
[13] B. Jourdain, T. Lelievre, and C. Le Bris. Existence of solution for a micro-macro model of 

polymeric fluid: the FENE model. J. Funct. Anal, 209(1):162-193, 2004. 
[14] B. Jourdain, T. Lelievre, C. Lc Bris, and F. Otto. Long-time asymptotics of amultiscalc model 

for polymeric fluid flows. Arch. Ration. Mech. Anal, 181:97-148, 2006. 
[15] D. J. Knezevic, and E. Siili. Spectral Galcrkin approximation of Fokkcr-Planck equations with 

unbounded drift. ESAIM: M2AN, 42(3):445-485, 2009. 
[16] A. Kufner. Weighted Sobolev Spaces. A Wiley-Interscience Publication. John Wiley & Sons Inc., 

New York, 1985. Translated from the Czech. 
[17] F. Lin, P. Zhang, and Z. Zhang. On the global existence of smooth solution to the 2-D FENE 

dumbbeU model. Comm. Math. Phys., 277(2):531-553, 2008. 
[18] F.-H. Lin, C. Liu, and P. Zhang. On a micro-macro model for polymeric fluids near equilibrium. 

Comm. Pure Appl Math., 60(6):838~866, 2007. 
[19] F. H. Lin and P. Zhang. The FENE dumbbell model near equilibrium. Acta Math. Sin. (Engl 

Ser.), 24(4):529-538, 2008. 
[20] P.-L. Lions and N. Masmoudi. Global existence of weak solutions to some micro-macro models. 

C. R. Math. Acad. Sci. Paris, 345(l):15-20, 2007. 
[21] C. Liu and H. Liu. Boundary conditions for the microscopic FENE models. SIAM J. Appl 

Math., 68(5):1304-1315, 2008. 
[22] H. Liu and J. Shin. The Cauchy-Dirichlet problem for the FENE model of polymeric flows. In 

Preparation. 

[23] T. Li and P. Zhang. Mathematical analysis of multi-scale models of complex fluids Commun. 
Math. Sci., 5:1-51, 2007. 

[24] N. Masmoudi. Well-posedness for the FENE dumbbell model of polymeric flows. Comm. Pure 

Appl Math., 61(12):1685-1714, 2008. 
[25] J. Necas. Sur une methode pour rcsoudre les equations aux derivees partielles du type elliptique, 

voisine de la variationnelle. Ann. Scuola Norm. Sup. Pisa (3), 16:305-326, 1962. 
[26] H. Ottinger. Stochastic Processes in Polymeric Liquids. Springer- Verlag, Berlin and New York, 

1996. 

[27] M. Renardy. An existence theorem for model equations resulting from kinetic theories of polymer 

solutions. SIAM J. Math. Anal, 22(2):313-327, 1991. 
[28] H. Zhang and P. Zhang. Local existence for the FENE-dumbbell model of polymeric fluids. 

Arch. Ration. Mech. Anal, 181(2):373-400, 2006. 

Iowa State University, Mathematics Department, Ames, IA 50011 
E-mail address: hliu@iastate.edu 

Iowa State University, Mathematics Department, Ames, IA 50011 
E-mail address: jmshin@iastate.edu 



